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Abstract 

Savannas are characterized by a discontinuous tree layer superimposed on a 
continuous layer of grass. Identifying the mechanisms that facilitate this tree- 
grass coexistence has remained a persistent challenge in ecology and is known 
as the "savanna problem" . In this work, we propose a model that combines a 
previous savanna model (Calabrese et al., 2010), which includes competitive 
interactions among trees and dispersal, with the Drossel-Schwabl forest fire 
model, therefore representing fire in a spatially explicit manner. The model 
is used to explore how the pattern of fire-spread, coupled with an explicit, 
fire-vulnerable tree life stage, affects tree density and spatial pattern. Tree 
density depends strongly on both fire frequency and tree-tree competition al- 
though the fire frequency, which induces indirect interactions between trees 
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and between trees and grass, appears to be the crucial factor controlling the 
tree-extinction transition in which the savanna becomes grassland. Depend- 
ing on parameters, adult trees may arrange in different regular or clumped 
patterns, the later of two different types (compact or open). Cluster-size 
distributions have fat tails but clean power-law behavior is only attained in 
specific cases. 

Keywords: Savanna, tree-tree competition, tree-grass equilibrium, 
individual based model, clustering, fire-spread model 

1. Introduction 



Savanna ecosystems are characterized by the robust coexistence of trees 
and grass. The mechanisms allowing for the persistence of both types of 
vegetation and governing the population dynamics and spatia l arrangement 



of sa vanna trees are poorly understood ( IScholes and Archer 



19971 : 



Bondl . 



20081 ). Of the many potential driving mechanisms investigated, local-scale 
i nteractions among trees have received increasing attention in re ce nt years 



Barot et al. 



2007 



1999; 



Wieeand et al. 



Calabrese et al. 



2006; 



Meyer et al. 



2008 



, 



2007alJbl;lScanlon et al. 



20101 ). Such tree-tree interactions can roughly be di- 



vided into two classes: facilitative and competitive. Facilitation among trees 
promotes tree clustering and may be mediated by a variety of mechanisms 
such as limited-range dis persal, improvemen t of local resour c e conditions, 



and protection from fire 



2005 



Scanlon et al 



20071 : 



Belsky et al 



1989 



Calabrese et al 



Hochberg et al 



2010) 



1994; 



Holdo . 



Alternatively, competition among trees for water, nutrients, and light may 
constrain tree density and favor tree-grass coexistence, as well as promot- 



ing spatial sepa r ation between trees (IBarot et al.l . Il999t iMeyer et al.l . 12008 ; 



Calabreseetal 



20101 ). 



There is evidence for both classes of interactions in the savanna litera- 
ture, sometimes coming from the same region. For example, several studies 
have 



1991 



bund evidence con s istent with competition in t 



Jeltsch et al. 



1999; 



Moustakas et al. 



2006 



r e Kalahari (Skarpe, 



2008 



Mever et al. 



while others have fo und evidence suggesting facilitation (jCaylor et al. 



2008) 



2003 



Scanlon et al. 



2007|). Indeed, one of the key difficulties in understanding the 



forces structuring savanna tree populations is that both classes of local-scale 
interactions often occur together and it is not obvious whether the net effect 
of loca l inter actions on tree population dynamics will be positive or negative 
( JBondl . 120081 ) . Further studies, both empirical and theoretical, are needed to 
better understand the interplay between these opposing forces. Specifically, 
studies that focus on a limited number of processes and their interactions 
should help illuminate the conditions under which positive or negative local 
interactions structure savanna tree populations. 

Mesic savannas that receive 400-800 mm of mean annual precipitation 
(MAP) are particularly interesting because there is evidence from such sys- 
tems t hat, in addition to local-scale interactions, fire plays an important 



role. (ISankaran et al. 



2005 ; 



Bucini and Hanan 



20071 ). Both of these fac- 



tors can act strongly on juvenile trees and can contribute to a demographic 
bottleneck through which juvenile trees must pass to recruit into the adult 
population. I n contrast to forest tre e species, savanna trees are often more 



fire resistant (JHoffmann et al. 



20031 ). thus savanna fires effectively burn the 



grass layer and the young trees included in it, leaving adult trees alive, af- 



fecting only tree recruitment and not adult survival ( JGignoux et al.l . 119971 ). 



Recent studi es highlighting the importance of tree competi ti on and/or fire on 



savannas are 



(2008) 



Higgins et a 



D'Odorico et al 



Calabrese et al 



(2000) 



Moustakas et al 



fl2006f ). lHananet al 



f 20 061) . JMoustakas et al 



fl2008f ). Mever et al.l ( 120081 ). or 



(120101 ). From their results we might expect a kind of tug of 
war between these forces, the outcome of which affects both the tree-grass 
balance of the savanna and the spatial arrangement of adult trees. 

The role of fire in mesic savannas is two-fold. On the one hand, it pro- 
vides an indirect way for grass to compete against trees: the higher recovery 
rates of grasses compared to juvenile trees make grass the dominant form of 
vegetation shortly after a fire has destroyed both. On the other hand, several 
studies have suggested that adult trees can pr otect vulnerable juveni l es from 



1994 



Holdo . 



fire, thus increasing their chances of survival (jHochberg et al. 
20051 ). but this protection effect has not been intensively studied. However, 
given the frequent occurrence of fires in many savannas, it seems likely that 
the protection effect may be one of the most common facilitative interac- 
tions among savanna trees, and the dominance of grass after fire could be as 
important as tree-tree competition in restricting the amount of tree-cover in 
the savanna. 
Recently 



Calabrese et al. 



(120101 ) studied the interaction between compe- 
tition and fire in a highly simplified savanna model. They showed that these 
two forces interact non-linearly with sometimes surprising c onsequences for 



Calabrese et al. 



tree p opulation density and spatial pattern. However, because 
(120 101 ) treated fire in a non-spatially explicit manner, only the negative im- 
pact on trees, and not the protection effect, was included and thus they 



could not fully tease apart how these contrasting local interactions function 
in combination. 

Here, we focus on a spatially explicit lattice model of savanna tree and 
grass population dynamics under the influence of competiti on and fire. The 
mode 



Calabrese et al 



is an extension of the semi-spatial model studied by 
(J2010|). Importantly, both competition and fire are spatially explicit pro- 
cesses in the new model. This allows us to study directly how adult trees 
influence the survival prob abilities of nearby juve niles. We treat competi- 



tion in the same way as in 



Calabrese et al. 



(120101 ) and fire is implemented 



in a similar manner as in the Drosse 



-Schwabl forest fire model from sta- 



tistical physics ( IDrossel and Schwab! 



19921 ). In contrast to adult trees in 



the Drossel-Schwabl model, grasses and juvenile trees are the flammable ob- 
jects in our case. We highlight the ranges of conditions under which local 
interactions result in net positive and net negative influences on juvenile 
tree recruitment, and we demonstrate how these local interactions affect the 
density and spatial structure of adult-tree populations. 



2. Spatially explicit fire models 



Bak and Chen! (119901 ) introduced a simple forest fire m odel to demon 



Drossel and Schwabl 



strate the emergence of scaling and fractal energy dissipation. 
( 119921 ) extended this model by introducing a lightning or sparking parame- 
ter /, and this is the forest fire model we have adapted to study fire spread 
in savannas. It is one of the best studied examples of non-conservative 



self-organized critica 



Drossel and Schwabl 



i tv (IBak and Chenl. 



1992 



Clar et al. 



1990 



1996 



Grassberger and Kantz 



1999; 



Schenk et al. 



1991 



20001). The 



forest fire model is a probabilistic cellular automaton defined on a 2-dimensional 
lattice of L 2 sites, initialized with a combination of burning trees and live 
trees, and updated at each time-step with the following four simple rules: (i) 
A burning tree becomes an empty site, (ii) A live tree becomes a burning 
tree if at least one of its nearest neighbors is burning. Some immunity can 
be introduced in this rule, so that a green tree becomes a burning tree with 



probability 1 — I (IClar et all 119961 ) . (iii) A new tree establishes at an empty 
site with probability r. (iv) Live trees in the lattice spontaneously (i.e., with- 
out the need of a burning neighbor) ignite with probability /. This model 
displays very rich behavior, and depending on the parameters / and r, it 
features spiral-like fronts, critical states and phase transitions. Furthermore, 
while the Drosel-Schwabl model is minimalistic, it produces burn patterns 



similar to those observed empirica 



wildfire models ( Zinck and Grimm 



l y, and is closely related to more detailed 



2009|). 



3. Savanna Fire Model (SFM) 

Our model is run in a square lattice with periodic boundary conditions. 
We use a lateral size of L = 200 sites, so that there are N = LxL = Ax 10 4 
lattice sites in the simulation domain. Each site represents a savanna square 



of 5 m eters on a side. In the previous savanna model (SM) of 



Calabrese et al 



(120 lOl ). each site in the lattice could be in one of two states: grass- or tree- 
occupied. The savanna fire model (SFM) introduced here combines the pre- 
vious SM and the above- described Drossel-Schwabl Forest Fire model, but 
with the flammable components being grass and juvenile trees. In this way, 
fire is included explicitly as a possible state in the dynamics. The SFM con- 



6 



siders three new states in addition to the two in the SM so that each site 
on the lattice can be in one of the following five states: Grass (G), Juvenile 
Tree (JT), Adult Tree (AT), Burning (B) and Ashes (A). 

We can distinguish two interaction neighborhoods for each lattice site: 
the near neighborhood consists of the eight sites sharing an edge or a corner 
with the central one (Moore neighborhood), and we assume this is the spatial 
scale at which direct competition among trees occurs. The far neighborhood 
consists of the sixteen additional sites surrounding the near ones and sharing 
edges or corners with them. They will be assumed to be the farthest sites to 
which seeds from a focal tree can arrive. 

We note that fire propagation occu rs over a much shor t er tim escale (the 



spread rate may be around 2 m/s, see ICheney and Gouldl (119951 )) than tree 



growth, reproduction, death, and other ecological processes. Thus we im- 
plement the burning process on top of the previous SM, but acting on a 
faster scale. Specifically, at each time step, time advances by At = 0.1 years, 
and the whole lattice is scanned in parallel to check for one of the following 
updates: 

1 . Growth: A random number is drawn for each site occupied by a juvenile 
tree so that with probability mAt it becomes an adult tree. Thus m~ l 
is the mean time for a juvenile tree to become adult. 

2. Reproduction and establishment: Each adult tree in the lattice sends, 
with probability bAt/24, a seed to each of the 24 sites within its near 
and far neighborhood . If the seed lands on a site in a state which is 
neither G nor A, then nothing happens (establishment fails). If instead 
a site occupied by grass or ashes is reached a juvenile tree is established. 



adult trees with probability P^ urv . This survival probability depends 
only on the competition exerted by neighboring adult trees: P^ u " 



3. Competition: A juvenile tree survives competition with neighboring 

dSi 

c 
. by neighboring adult trees: 1 c 

e~ SSl , where 8 is the competition parameter and Si is the number of 

adult trees in the near neighborhood. 

4. Death: A random number is drawn for each site occupied by an adult 
tree, so that with probability aAt this tree dies. Thus a^ 1 is the 
average adult-tree lifespan. 

5. Recovery: At each time step, each ash site may recover into grass with 
probability rAt, so that r" 1 is the mean recovery time of grass from 
ashes. Note that this forces a delay between successive fire fronts, thus 
preventing the lattice from continually burning. 

6. Spontaneous burning: There is a "lightning parameter", /, so that fire 
appears spontaneously on the lattice at this rate, affecting grass and 
juvenile trees. More explicitly, lattice sites occupied by G and JT are 
checked so that with probability fAt/N they become burning sites. 

7. Fire propagation and extinction: After updating with all the above 
processes, a new pass through the lattice is done, so that if some fire 
has been introduced in the previous step, fire propagation is simulated 
until the fire burns out. As previously mentioned, we assume that 
this process is fast and occurs on a much shorter timescale than the 
At = 0.1 years introduced above. It is implemented in the following 
way: 

a) Each G and JT site is checked and if at least one site in its near 
neighborhood is in the B state, the site also burns with probability 



1 — I, where I & [0, 1] is an immunity parameter. This models fire 
propagation on grass and juvenile trees. Note that, since adult 
trees do not burn, fire has a lower chance of reaching JT (and 
G) sites which are surrounded by some adult trees. The inclu- 
sion of fire in a spatially explicit manner therefore implements the 
protection effect from adult trees. 

b) End of burning: All sites that were burning before entering the 
previous step (a) are set to ashes. 

Processes a) and b) are repeated until no burning sites remain. Time 
then advances At units and the algorithm repeats again from step (1) 
on the updated lattice. 

Tabic 1: Parameters 

Parameters Units 

a adult-tree mortality rate 0.01 year~ l 

b seed dispersal rate 0.05 year' 1 

5 tree competition coefficient vary 

/ lightning parameter 0.33 year -1 

m juvenile tree maturation rate 0.1 year -1 

r grass recovery rate 4 year -1 

I fire immunity 0.3 



The core parameters used (see Tabled]) are based on those of 



Calabrese et al 



( 120101 ). but with a few modifications. The adult-tr ee death rate , p, an d the 



lightning parameter, /, were changed according to 



Hanan et al. 



(120081 pag 



852) and iGignoux et al.l (119971 . pag. 557), respectively. In m esic savannas 
fire fre quency is abou t once per year t o onc e every three years (IHanan et al 



(12008 



pag. 852) and iGignoux et al.l (J1997I . pag. 557) ) so we set f = 0.33 



Hochberg et al. 



(1994. 



The juvenile tree growth rate, m, was determined from 
pg. 219), so that in the absence of fire, a juvenile tree takes on average 5 
years to reach the adult state. These parameters will be used through the 
paper unless explicitly stated. 

The spatially-limited dispersal of seeds from an adult to neighboring sites, 
one of two key facilitative processes in the model, occurs at the spatial scale 
of the first and second neighborhood, whereas the main negative interac- 
tion, adults competing with and inhibiting the development of juvenile trees, 
occurs only at the scale of the first (near) neighborhood. 

This is the opposite situation as the one believed to occur in extremely 
arid ecosystems, namely short distance facilitation by local improvement 
of water infiltration, and long range competition among plants mediated 
by long superficial roots. In this last case, vegetation p atterns are ex- 



pected to display a rather regular tree or patch spacing ( IRietkerk et al 



2004 



Rietkerk and van de Koppel 



20081 ). Our situation is more appropriate 



for mesic savannas, and would tend to promote tree clustering. But the oc- 
currence of fire may alter the nature of the interactions in a variety of ways, 
which we investigate in the following. 

4. The tree-grass balance and tree extinction 

As expected, stronger tree-tree competition shifts the tree-grass balance 
in favor of grasses (Figure HJ left). The model was run for 5000 years to ensure 



10 



an asymptotic state was reached in which we performed the measurements 
described in Figure [TJ Simulations were performed to determine under which 
conditions a transition fro m savannas to grasslan d occurs. 



Similarly to the results 



Calabreseetal 



( J2010l ) obtained for their fire pa- 



rameter a, the lightning frequency / turns out to be the parameter with 
the strongest influence on the savanna-grassland transition in the SFM: The 
right part of Figure [I] shows a phase t ransition from savanna to grassland 



Calabrese et al. 



(boioh ). 



driven by increasing / (c.f. Figure 2 in 

Frequent fires prevent juvenile trees from recruiting into the adult pop- 
ulation, and if this inhibition is strong enough, it would eventually result in 
tree extinction. This mechanism and the subsequent tree extinction it causes 
was implicitly c o ntain ed in the definition of the fire parameter in the SF of 



Calabrese et al. 



(I2010f l. Here, the mechanism appears as a consequence of 



the explicit presence of fire. 

5. Positive and negative effects of surrounding adult trees on ju- 
venile's: Protection vs. competition 

In addition to affecting the tree-grass ratio, fire also introduces the posi- 
tive effect of juvenile tree protection by surrounding adult trees. To analyze 
this effect in detail we ran simulations in which only the fire propagation 
process (step 7 in the above algorithm) occurs. The lattice is initialized fully 
with grass except for one unique site occupied by a juvenile tree, and a num- 
ber of adult trees, from 1 to 8, occupying random positions in the Moore 
neighborhood of the JT. Given this initial condition, one sparking is allowed 
so that a lattice site chosen randomly among the G sites burns and fire begins 



11 



to propagate. Step 7 in the SFM algorithm is repeated until fire disappears. 

Juvenile trees sufficiently protected by adult trees will not burn. An 
example is seen in Fig. |2j where the pass of a fire front does not affect a 
juvenile protected by five adult trees. 

We quantified this effect by repeating the burning protocol 1000 times 
for each number of AT neighbors, from 1 to 8. The sparking site and the 
position of the surrounding neighbors is changed randomly in each of these 
realizations. The resulting survival probability Pp urv (Si) is shown in Fig. 
|3l The protection provided by an increasing number of ATs is clearly seen 
when the immunity parameter is not too small. For very small /, protection 
is only effective when the juvenile is completely surrounded by adults, i.e. 
S 1 = 8. 

To better quantify the impact of the protection effect on juvenile survival 
and recruitment, we now estimate how the number of adult trees Si in the 
near neighborhood of a site affects the recruitment probability -Pr(Si) defined 
as the probability that a grass site becomes successfully colonized by a tree 
seed during a given time-step and the resulting JT survives successive fires 
and becomes an adult. This probability is a product of several factors. 

First, the grass site should receive in that time step a seed from the 
adult trees in the near or in the far neighborhood (the numbers of adult 
trees there are Si G (0,8) and S2 G (0,16), respectively). This is given by 
P s (Si, S 2 ) = 1 — (1 — bAt/24) Sl+S2 . Then, the seed establishes as a juvenile 
tree and must survive competition during successive time steps, which is 
given by the factor P^ urv = exp (— 8 Si). Since m _1 is the average growth 
time from juvenile to adult tree, (mAt) -1 time steps occur during growth, 



12 



and exp (— -^i) is the total survival factor to adulthood under competition. 

Finally, the growing JT should resist the first and successive fires occur- 
ring during its growing time m~ l . The probability of surviving a single fire 
is the function Pp urv (Si) numerically calculated and shown in Fig. [3] for the 
case in which the focal site is surrounded by S\ trees in the near neighbor- 
hood (i.e., S2 = 0). An estimation of the probability surviving successive 
fires, which neglects any correlations arising from successive fire fronts and 
from the lattice configuration beyond the immediate neighborhoods, would 
be (Pp urv (Si)) m , where f/m is the expected number of fires suffered by the 
JT during its growing time m~ l . The probability P s {Si, S2) depends on the 
number of AT both in the near and in the far neighborhood. For consistency 
with the calculation of Pp v (Si) we will take S2 = 0. This (as in the case 
of Pp v (Si)) will underestimate the probability of establishment, survival 
and recruitment, as the trees in the far neighborhood do not compete with 
the central one. In this way we will obtain an estimation of the recruitment 
probability Pr(Si) that is smaller than the exact one. Thus, if this function 
shows positive effects of surrounding adult trees, the exact result must be 
larger, since our approximation is obtained in a worst case situation. 

Summarizing all the factors above, with 5*2 = 0, our estimation of the 
recruitment probability of a grass site surrounded by Si adult trees is 



Pr(Si 



1 



6At x Sl 



e-^(p|«™(S 1 ))- (1) 



24 

This is plotted in Fig. HJ and reveals both the positive and the negative 
effects of the presence of neighboring trees (but remember that the positive 
effects are underestimated). For medium values of the competition parameter 

13 



and above four neighboring adult trees the positive protective effect of fire 
(in combination with local dispersal) overcomes the negative effect of direct 
competition (see Figure H^a) and (c)), but for high values of competition the 
negative effect predominates (see Figure ID^b)). For frequent fire, however, 
the protection effect is no longer effective. 

6. Clustering patterns 

6.1. Tree spatial pattern under different fire scenarios 



We c haracterize spatia 



function ( IDieckmann et al. 



patte rns of adult trees by the pair correlation 



2m,g(l): 



9(1) 



Paa{1) 
(Pa) 2 



(2) 



where paa is the proportion of pairs of adult trees at a distance I (with 
respect to the total number of pairs of sites at that distance) and the denom- 
inator is the expected value of this proportion under a random distribution 
with the density of the adult trees pA- At large distances g(l) is expected to 
approach 1, as correlations indicating a departure from random distribution 
would decay. For short distances, q(l) characterizes how the trees are packed 



together (see 



Dieckmann et al 



2001 



, chap. 14), values higher than 1 indi- 



cating a proportion of pairs at that distance greater than in the random case 
(clustering), and a smaller proportion indicated by values of g smaller than 
1 (revealing a more regular spacing). We will not use the Euclidean distance 
for I but instead we will measure / in number of cell layers so that g(l) and 
g(2) will denote the pair correlation function for the first and for the second 
Moore neighborhood, respectively. 
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Comparison of Figure \5\ with the results of ICalabrese et al.l (120101 ) shows 
that all the patterns found in the SM are also present here. Some features 
of the patterns can be understood from the fact that there is direct compe- 
tition only between nearest neighbors, whereas the facilitation effect of local 
seed dispersal reaches first and second neighbors. In consequence, all these 
patterns have an enhanced probability of ATs having other ATs as second 
ne ighbors (far Moore nei ghborhood), as seen by the high value of g(2). As 



m 



Calabrese et al. 



(120101 ). two types of configurations are distinguished by 
having a value of g(l) smaller or larger than 1, i.e. smaller or larger propor- 
tion of ATs in the near neighborhood than the one expected from a random 
distribution. The balance between positive and negative tree-tree interac- 
tion effects determines these values. The case g(l) < 1 is a regular case in 
which trees appear more regularly spaced than in the random case. The case 
g(l) > 1 is a clumped state, in which, although the density of near-neighbor 
pairs is still smaller than the one of far-neighbor pairs, it is larger than in 
the random case. The transition between the two states was governed by 
the parameter a in the SM, which controls the probability of surviving fire. 
Here, this transition is determined by the explicit fire parameter /. 

In the clumped patterns just described, further illustrated by Fig. E^a) 
and (b), the clusters are open in the sense that there are more neighbors in 
the far neighborhood than in the near neighborhood. This a clear effect of 
the competition existing in the near neighborhood, and was the only clumped 
state present in the previous SM. The novelty here is that, in addition, there 
is a second type of clustered state not present in the SM. A clumped state 
made of closed clusters is illustrated by figure [6(c) and (d). The clusters are 
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closed in the sense that there are more AT neighbors in the near neighborhood 
than in the far neighborhood. Thus, the positive effect of fire protection (and 
local dispersal) has completely overcome the competition effect occurring in 
the near neighborhood. The transition from one type of pattern to the other 
occurs when changing the competition or the lightning parameters, S and /, 
as shown in Figure [71 



6.2. Cluster-size distributions 

A cluster is a group of neighboring sites occupied by the same type of 
vegetation (e.g. adult trees). The distribution of cluster size s is a powerful 



indicator of the different mechanism s occurring in ecosystems (jPascual et al. 



2002 



Pascual and Guichard 



20051 ) . A dult tree clust e r-size distributions in 



the Kalahari have been investigated by 



Scanlon et al. 



(120071 ). finding that in 



most cases a power- 
of uniform quality) . 



aw fit can describe the data (although the fit was not 



Scanlon et al. 



( 120071 ) showed that resource constraints, 
together with positive local interactions of the type identified in the previous 
section, could generate cluster-size distributions similar to the observed ones. 
Figure [8] shows complementary cumulative distributions of adult-tree clus- 
ter sizes from our model, where the Moore neighborhood has been used to 
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define clusters. Though the distributions have fat tails (see Fig. [S]), a sin- 
gle power law does not provide a good description in the realistic range of 
parameters considered above. Also, the plateau at large sizes in the small- 
/ curves of Fig. [8] indicate the presence of clusters much larger than the 
rest. By artificially changing parameters to other ranges, one can find situ- 
ations in which the cluster-size distribution follows a relatively good power 
law. This happens, for example, for a = (3 = 1, S = 0.01, and / « 0.9. 
Inspection of the tree distributions above and below this /-value indicates 
that a percolation transition occurs precisely at that point: there is a giant 
AT cluster spanning the whole area for smaller / values, and disconnected 
tree patches for higher values. Power-law cluster distributions are observed 
then close to this percolatio n transition, as in the mechanism discussed in 



Pascual and Guichardl ( 120051 ). although in a narrower parameter range than 
suggested there. This transition is not attained within the parameter ranges 
considered before in this paper: tree cover in Fig. [S]is just of 0.4 for / fh 0.9, 
and can not be increased much more (see Fig. [TJ right panel), which makes 
it difficult to attain percolation through the whole lattice because of the ab- 
sence of very large clusters. By artificially changing parameters to obtain 
larger tree densities, percolation becomes easier. In such situations, we ob- 
serve more robust power-law behavior (not shown), but the system is then 
closer to a forest than to a realistic savanna. We do not find systematic 
correlation between the small-scale character of the tree patterns (regular, 
clumped, open, closed, ...) and the type of cluster-size distributions, despite 
the fact that one could expect that p ositive short-range correlations would 



favor power-laws (jScanlon et al. 



20071). 
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7. Summary 

We have introduced a model for savanna structure which includes, in ad- 
dition to fundamental ecological interactions including competition, the ef- 
fects of spatially explicit fire spread. Fire introduces some effective tree-grass 
and tree-tree interactions which are important in shaping tree demography 
and spatial pattern. First, the presence of fire improves competitiveness of 
grass because of its faster post-fire recovery. Second, adult trees may protect 
nearby juveniles from fire. This results in a positive tree-tree interaction 
which can, in some circumstances, overcome the effects of tree-tree compe- 
tition for resources. A variety of tree spatial distributions are observed as 
a result of these direct and indirect interactions, which we have character- 
ized by the pair correlation function and the cluster-size distribution. As 
the short-range positive interactions gain importance relative to the negative 
ones, a succession of regular to clumped states is observed. Clumped states 
can have "open" clust ers, like the ones present in the previous SF model 
( jCalabrese et al.l . 120101 ) , but also "closed" clusters for the cases with stronger 
positive interactions. Adult-tree cluster-size distributions are of power-law 
type in some cases because of the proximity to a percolation transition, but 
for much of the realistic parameter range tree cover is small and far from 
percolating. The tails of the distributions, although fat, seem to decay faster 
than power laws, as seen in fact to occur in several of the sites reported by 
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(120071 ). 
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Figure 1: Density of adult trees (i.e., the number of adult trees divided by the total number 

of lattice sites) versus competition (J, left graph) and versus lightning (/, right graph). 

Average over 500 snapshots in the long-time asymptotic state. Parameters are as in Table 

[T] Transition from the coexistence state to grassland is driven by increasing S and/or /. 
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Figure 2: Protection effect: Selected snapshots from an example simulation with 5 adult 
trees (blue) surrounding a juvenile (green). Immunity parameter / = 0.3. Time runs 
from left to right and then from the upper to the lower row. A fire front (red) advances 
downwards, converting grass (yellow) into ashes (white), but the juvenile survives. Only 
a 10 x 10 area of the whole 200 x 200 lattice is shown. 
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Figure 3: The protection effect: The probability Pp urv (Si) of a juvenile surviving one 
fire as a function of the number of surrounding adult trees in its first neighborhood. This 
probability has been obtained from 1000 realizations of the process in which fire is initiated 
at one grass site, as described in the text, using immunity / = 0.3. 
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Figure 4: Estimation of the recruiting probability Pr(Si), as a function of the number of 
adult trees Si in the near neighborhood, from Eq. ([1]), showing the positive and negative 
effects of these neighbors, (a) and (b) / = 0.3, / = 0.33 year" 1 (triennial fire), (c) / = 0.3, 
/ = 0.2 year -1 (pentannual fire). 
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Figure 5: Patterns in the SFM. Parameters as in Table Q] and S — 0.01. Regular case: 
/ = 0.19. Clumped case: / = 0.31. The central panel shows an intermediate state 
(/ = 0.223) in which <?(1) = 1, which indicates the same number of AT near pairs as in a 
random case. 
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Figure 6: Savanna configuration in the clumped state at parameters as in Table [TJ Panels 

(a) and (c) display only a 30 x 30 portion of the full 200 x 200 lattice, (a) S = 0.01, 
/ = 0.31: an example of open cluster of ATs, the typical configuration at these parameters. 

(b) Pair correlation function, similar to the one in the SM in the clumped state, (c) 
5 = 0.001, / = 0.45: an example of closed cluster of ATs, the typical configuration at 
these parameters, (d) Pair correlation function, which is different to the one in the SM 
in the clumped state because the maximum of g(l) occurs at I = 1, i.e. in the near 
neighborhood. 
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Figure 7: (a) Values of the pair correlation function for near neighbor pairs g(l) and for far 
neighbors g (2), and values of adult-tree density, both as a function of the tree-competition 
coefficient 8. Closed clusters occur for 8 < 5.2 • 10~ 5 and open ones for 8 > 5.2 • 10 -5 . (b) 
8 = 0.001, Pair correlation functions showing the transition between open (5(1) < 5(2)) 
and closed (g(l) > g(2)) clusters by increasing /. Other parameters in both panels as in 
Table [U 
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Figure 8: Left: Adult-tree cluster-size distributions, represented by means of the comple- 
mentary cumulative distributions. Parameters as in Table[T]but 8— 0.001. Right: Density 
of adult trees (i.e., the number of adult trees divided by the total number of lattice sites) 
versus lightning /. 
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